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ABSTRACT 

Saturn's C ring contains multiple spiral patterns that appear to be density waves driven by periodic 
gravitational perturbations. In other parts of Saturn's rings, such waves are generated by Lindblad 
resonances with Saturn's various moons, but most of the wave-like C-ring features are not situated near 
any strong resonance with any known moon. Using stellar occultation data obtained by the Visual 
and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft, we investigate the origin 
of six unidentified C-ring waves located between 80,900 and 87,200 km from Saturn's center. By 
measuring differences in the waves' phases among the different occultations, we are able to determine 
both the number of arms in each spiral pattern and the speeds at which these patterns rotate around 
the planet. We find that all six of these waves have between 2 and 4 arms and pattern speeds between 
1660°/day and 1861°/day. These speeds are too large to be attributed to any satellite resonance. 
Instead they are comparable to the predicted pattern speeds of waves generated by low-order normal- 
mode oscillations within the planet [Marley & Porco 1993, Icarus 106, 508]. The precise pattern 
speeds associated with these waves should therefore provide strong constraints on Saturn's internal 
structure. Furthermore, we identify multiple waves with the same number of arms and very similar 
pattern speeds, indicating that multiple m = 3 and to = 2 sectoral (/ = ra) modes may exist within 
the planet. 

Keywords: planets and satellites: interiors — planets and satellites: rings 



1. INTRODUCTION 

Spiral waves are patterns in Saturn's rings produced 
by periodic gravitational perturbations on the ring ma- 
terial. Many of these waves can be attributed to mean- 
motion resonances with Saturn's various moons. How- 
ever, a number of spiral waves in the C ring do not fall 
near any known resonance with any moon. These waves 
were first n oticed in the Voyager radio occultation data 
(jRosen et a l. 1991), but have also been observed in stel- 
lar occultations by the ultraviolet spectrometer onboard 
the C assini spacecraft ([Colwell et al.ll2009bl : iBaillie et al.l 
120111) . However, these previous studies were unable to 
identify the source of the perturbations responsible for 
generating these waves because they could not determine 
either the number of arms in the spirals or how fast the 
relevant patterns were rotating around the planet. Using 
occultation data from the Visual and Infrared Mapping 
Spectrometer (VIMS), we have now been able to deter- 
mine both these quantities for six of these unidentified 
waves. These six waves have the right pattern speeds 
and symmetry properties to be produced by low-order 
normal mode oscil lations within Saturn, a s predicted by 
iMarlevI ()1990[ ) and lMarlev fc Pored (|1993[ ). These waves 
should therefore provide valuable new constraints on the 
planet's interior structure 

Prior to describing our analysis of these waves, we first 
provide some background information about the relevant 
spiral waves in Section [2] This includes a brief summary 
of the theory behind spiral waves (Section 12. 1|) . a sum- 



mary of the wave-like features in the C ring fSection l2.2p 
and a description of the six waves that will be investi- 
gated here (Section l2.3p . Section [3] then describes how 
we can determine the number of arms and pattern speeds 
of these waves. Section [3T] describes the VIMS occulta- 
tion data that will be used in this analysis. Section 13.21 
illustrates how comparing multiple occultations can con- 
strain the waves' symmetry properties using a particu- 
larly informative set of occultations. Sections [331 and [3^ 
describe the wavelet-based techniques we employ to as- 
certain the number of arms and pattern speeds of these 
waves. Section |3] summarizes the results of our calcula- 
tions. Finally, Section [5] discusses the implications of our 
findings. 

2. BACKGROUND 

2.1. The theory of spiral waves 

Spiral patterns in Saturn's rings are generated by vari- 
ous periodic perturbations on the ring material, and dif- 
ferent types of perturbations produce different types of 
spiral waves. For example, periodic vertical forces gen- 
erate warped structures known as bending waves, while 
periodic radial or azimuthal forces produce variations in 
the rings' surface density known as density waves. The 
structure and dynamics of these spiral patterns are rela- 
tively well understood and good theoretical overviews of 
these phenomena are available (e.g. Shu 1984). Hence, 
we will only briefiy review the aspects of spiral waves 
that are most relevant to this analysis. Note that the 



waves considered below all appear to be density waves, 
so we will focus on the dynamics of those structures here. 

A given spiral density wave consists of an integer num- 
ber of arms that become more tightly wrapped with in- 
creasing distance from the radius of the exact resonance. 
This entire pattern rotates around the planet at a sin- 
gle pattern speed. The waves in Saturn's rings are so 
tightly wrapped that the opacity variations appear to 
be almost purely radial in both images and occultation 
profiles. Hence the waves appear as periodic variations 
in ring brightness or opacity whose wavelength changes 
systematically with radius (see Figure [1] below for some 
examples). However, because a density wave is actually 
a rotating spiral pattern, the locations of the peaks and 
troughs in a given profile will depend on both the ob- 
served ring longitude and the orientation of the pattern 
during the measurement. 

The organized motions responsible for a spiral density 
waves are most efficiently generated by Lindblad reso- 
nances. At these resonances, the ring-particles' radial 
epicyclic frequency k is an integer multiple of the differ- 
ence between the angular frequency of the perturbing po- 
tential rip and the ring-particles' mean motion n. Hence, 
if there is a density wave driven at a given radial loca- 
tion in the rings, the most likely resonant perturbation 
frequencies will satisfy the following relationship: 



lowing asymptotic expression 



m[n — Q,p) = K 
where m is any non-zero integer (i.e., m 



(1) 



.-3, -2, -f, 
1, 2, 3...). Rewriting the resonant condition in terms of 
the local apsidal precession rate, n7 = n — k, we have the 
familiar expression for a first-order Lindblad resonance: 



{m — l)n + 



77lflr. 



(2) 



Since w << n, resonances with ttt, > (known as inner 
Lindblad resonances or ILRs) have ftp « (m — l)n/m, 
while those with to < (known as outer Lindblad reso- 
nances or OLRs) have fip ss {\m\ + l)n/|?Ti|. 

In a differentially-rotating, self-gravitating disk, the 
periodic perturbations at such a resonance give rise to 
a trailing spiral density wave that propagates away from 
the location of the exact resonance towards the location 
in the rings where n = Qp. Thus, for a Keplerian disk 
like Saturn's rings (and assuming J7p > 0), density waves 
will propagate outwards from an inner Lindblad reso- 
nance and inwards from an outer Lindblad resonance. 
The pattern of surface density variations generated by 
such a wave gives rise to variations in the local optical 
depth with radius r, inertial longitude A and time t. For 
waves of small amplitude these variations may be written 
as: 

r (r. A, i) ~ To -h Ar (r) cos 4>{r,X,t), (3) 

where the wave's phase (fi can be decomposed into a part 
that depends only on the observed longitude and time, 
and another that depends only on radius: 



0(r,A,t) ~ \m\{X-npt) + (j>r{r). 



(4) 



Hence \m\ gives the number of arms in the wave pattern, 
while flp is the angular rate at which it rotates around 
the planet. The wave's pattern speed therefore equals 
the angular frequency of the external perturbing force. 

At sufficiently large distances from the resonance, the 
radius-dependent part of the phase is given by the fol- 



(j)r{r) 



3(m-l) + j4(^ 
2 \rL 



Msjr - tlY 



where Ms is the mass of Saturn, J2 is a measure of the 
planet's oblateness, rs = 60, 330 km, r^ is the resonant 
radius where Eq. ([T} is satisfied, ao is the undisturbed 
surface mass density of the ring and (po is a constant. 
Note that for an outward-propagating wave from an ILR 
both m — 1 and r — r^ are positive, while for an inward- 
propagating wave from an OLR m — 1 and r — r^ are 
both negative. Hence in both cases d(f>/dr > and a 
line of constant phase will have dr/dX < 0, correspond- 
ing to a trailing wave pattern. Also, the wave's radial 
wavenumber k derived from the above expression: 



fc(r) 



dr 



3(to - 1) + J2 



21 



Msjr - tl) 



(6) 
increases linearly with distance from r^. 

So long as the background opacity of the ring Tq does 
not vary too much with radius, the phase parameter (f> 
can be computed directly from the local opacity varia- 
tions {e.g., opacity maxima occur where </) ~ and opac- 
ity minima occur where (j) — 180°). Thus 0(r, A) can be 
calculated from any given radial profile across the ring. 
Indeed, for identifiable waves the trends in (f) with radius 
can be used to est imate parameters like the ring's sur- 
face mass density (jTiscareno et al.|[2007[) . However, for 
our efforts to determine the pattern speeds and m-values 
of unidentified waves, it is more useful to consider the 
difference in phase parameters between occultations ob- 
served at different times and longitudes. For these phase 
differences S(j), the radius-dependent term (j^r should can- 
cel out, leaving terms that only depend on the known 
differences in the observation times {6t) and observed 
longitudes (SX), and the unknown parameters \m\ and 

i Lp. 



\m\{dX-np6t). 



(7) 



We can therefore compare the observed phase differences 
5(f> with those predicted for various combinations of m 
and f2p, and thereby determine which pattern speed and 
m-number is most consistent with the observed phase 
differences. 

2.2. Spiral waves in the C ring 
Surveys of Saturn' s C-ri ng by iRosen et al.l ()199lD 



gg 

iColwellet ail ()2009bO and iBailhe et al.l ()2011[ )'" have 



found a total of 27 wave-like structures in the C-ring. 
Only five of these features have been identified with 
specific satellite resonances. The most prominent fea- 
ture is a bending wave generated by an unusu al Ti- 
tan -1:0 nodal res onance (jRosen fc Lissaueil 119880 . and 
IRosen et al.l (|199lD identified one density wave produced 
by the Mimas 4:1 ILR. More recently, the density wave 
produced by the Atlas 2:1 ILR has been detected, and 
there are hints of density waves at the lo cations of the 
Mimas 6:2 and Pandora 4:1 resonances (jColwell et al.l 
I2009bl: iBaillie et ai1l201lD . More than half of the remain- 
ing, unidentified waves appear to be inward-propagating 
{i.e., their radial wavenumber increases inwards). Such 
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behavior is inconsistent with density waves generated by 
resonances with Saturn's moons. Since all the moons or- 
bit outside the rings, the perturbation frequency associ- 
ated with any moon's gravitational perturbations will be 
smaller than the local orbital frequency anywhere in the 
C ring. Hence Lindblad resonances in the C ring with any 
moon should be ILRs, which should generate outward- 
propagating waves. The inward-propagating waves must 
therefore be produced by some other mechanism. 

One possible explanation for these inward-propagating 
waves is that they are bending waves instead of den- 
sity waves. Unlike density waves, bending waves typi- 
cally pr opagate aw ay from the location in the rings where 
flp = n (|Shulll98^ . Hence vertical resonances where the 
perturbation frequency is less than the local orbital fre- 
quency (the so-called inner Vertical resonances or IVRs), 
like those generated by several of Saturn's moons, should 
propagate inwards. Furthermore, the very low ring open- 
ing angle in the Voyager radio experiment {B = 5.6°) 
meant that the vertical warps in the bending wave could 
lead to significant optical depth variations along the ob- 
served line of sight. However, most of these inward- 
propagating waves are also clearly visible in stellar oc- 
cultations observed by Cassini at high ring opening an- 
gles, where optical depth variations from bending waves 
would be very subtleQ Furthermore, there are no strong 
candidate IVRs due to known (or even hypothesized) 
satellites at the desired locations. Thus this is not the 
currently- favored interpretation for most of these inward- 
propagating features. 

An alternative explanation of these waves is that they 
are not generated by resonances with Saturn's moons, 
but instead are produced by resonan ces with normal - 
mode oscillations within Saturn itself. iStevensonI (|1982[ ) 
first suggested that oscillations within the planet could 
produce identifiabl e str uctures in the r ings. Later, 
iMarlevI (|T990l HoMI ) and lMarlev fc Pored (fl99l demon- 
strated that acoustic modes in the planet's interior could 
indeed give rise to gravitational perturbations with pat- 
tern speeds that are sufficiently fast to generate OLRs 
and inward-propagating waves in the C ring. The lat- 
ter authors also found that the predicted pattern speeds 
of low-degree /-modes were consistent with some of the 
observed wave locations, but uncertainties in Saturn's in- 
terior structure and in the theory of normal modes meant 
that the resonance radii could only be predicted to within 
- 500 km. 

2.3. The six waves examined in this study 

This analysis will focus on six of the unidentified, 
inward-propagating waves. As summarized in Table (TJ 
these waves have been designated in var ious ways by 
different authors. Here we will use the iColwell et al.l 
(j2009b( ) nomenclature, which identifies each wave with 
a number giving the radial location of the wave in thou- 
sands of kilometers. Figure [T] shows profiles of all six 
waves derived from the highest resolution VIMS occul- 
tation to date (see below for how these data were pro- 
cessed) . Each wave appears as a periodic variation in the 



TABLE 1 

Summary of waves investigated here, including 
designations by various authors. 



Wave 


Rosen et al. 


Colwell et al. 


Baillie et al. 


Location" 


(1991)'' 


(2009b) 


(2011) 


80988 km 


e 


W80.98 


13 


82010 km 




W82.00 


15 


82061 km 


f 


W82.06 


16 


82209 km 


g 


W82.21 


17 


84644 km 


i 


W84.64 


19 


87189 km 




W87.19 


32 



^ One exception is the wave designated "j" in lRosen et al.l 11199 II ), 
which is in fact invisible in all VIMS occultation profiles obtained to 
date, but is seen in one or more UVIS stellar occultations observed 
at very low opening angles (Baillie et al. 2011 ). 



" Inferred resonance location obtained by fi tting Eg. [6] to UVIS 
occultation data, see Table 7 of Baillie et al] l l 201111 . 
^ See Figures 7 and 8 in lRosen et al. (_1991j). 



ring's opacity, with a wavelength that varies with dis- 
tance from the planet. More specifically, all these waves 
have wavelengths which decrease inwards, which suggests 
that they are either bending waves driven by IVRs or 
density waves driven by OLRs. Since the visibility of 
these waves does not appear to vary with spacecraft ele- 
vation angle, it is unlikely that any of these features are 
bending waves. Indeed, their pattern speeds turn out to 
be consistent with those expected for OLRs. Hence for 
the remainder of this discussion we will anticipate our fi- 
nal result by referring to these features as density waves. 

All six waves are located in the central part of the C 
ring, between radii of 80900 and 87200 km. This is a re- 
gion of gently undulating structure, with an average nor- 
mal optical depth {i.e., the optical depth the ring would 
have if the line of sight was exactly perpendicular to the 
ringplane) of about 0.10. All but the outermost wave are 
located between the 'plateau' features designated P4 and 
P5 by Colwell et al. (2009b). While four of the waves faU 
in otherwise unexceptional locations, two fall either on 
top of (W87.19) or immediately adjacent to (W82.06) 
local maxima in optical depth. 

The strongest known eccentric resonances in this region 
are the two 4th-order Enceladus 5:1 resonances at 82538 
and 82477 km, the 3rd-order Janus and Epimetheus 5:2 
resonances at ^ 82780 km and ^ 82950 km, the Pan 2:1 
ILR at 85105 km, and the Atlas 2:1 ILR at 87647 km. 
Only the last of these is known to be associated with 
an observable wave. None of these resonances is located 
within 100 km of any of the unidentified waves. There 
are two very weak vertical resonances in the vicinity of 
our target waves: the Enceladus 5:l(e^«) resonance at 
80967 km (located 21 km interior to the W80.98 wave), 
and the Janus 5:2(e«) resonance at 82096-82107 km (^^40 
km exterior to the W82.06 wave). However, as mentioned 
above, the visibility of these waves at large ring opening 
angles argues against these features being bending waves, 
and furthermore many other comparable (or stronger) 
vertical resonances do not produce visible bending waves 
in the rings. Hence we consider these rough alignments 
to be coincidental, and conclude that none of these six 
waves can be attributed to a resonance with any of Sat- 
urn's moons. 

These six waves have the strongest opacity variations 
and the longest wavelengths of any of the unidentified C- 
ring waves. We focus exclusively on these waves because 
they provide the best opportunities for determining un- 
ambiguous pattern speeds. The strong opacity variations 
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Fig. 1. — VIMS occultation profiles of the six waves examined in this analysis. Each panel shows the ring's normal optical depth versus 
ring radius, which is measured in kilometers from the inferred resonance locations r^. (For the innermost five waves, the r^ value comes 
from Baillie et al. 2011, while for W87.19 the resonance position has been adjusted to match the best-fit pattern speed of this wave, 
see Seciton |4]) The specific profiles shown here come from an occultation by the rings of the star R Cassiopaea, which provides our 
highest-resolution profiles of these waves. The raw data numbers were converted to transmission estimates by normalizing the stellar signal 
to unity in the middle of the Maxwell Gap (87,375-87,425 km), and then translated to normal optical depth values using the standard 
formula, assuming the elevation angle of the star is 56.04° above the ringplane. 



mean that we can clearly detect the maxima and min- 
ima, which makes the relevant phase parameters easier to 
determine. Furthermore, the relatively long wavelengths 
of these waves should minimize our sensitivity to small 
errors in the occultation geometry. Errors in Cassini's 
trajectory reconstruction can cause the occultation pro- 
files to shift slightly in radius, which can in turn create 
problems when comparing data from different occulta- 
tions. For example, if one profile was shifted relative to 
the other by half a wavelength, then the phase difference 
between the waves might appear to be zero when in re- 
ality it is 180°. However, this is unlikely to occur with 
any of these six waves, because they all have maximum 
wavelengths exceeding 2 kilometers, which is consider- 
ably larger than the sub-kilometer uncertainties in our 
reconstructions of occultation geometry (see below) . 

3. METHODS 

3.1. VIMS Observations 

The VIMS inst rument is described in detail in 
iBrown et all ()2004[ ). While VIMS is typically used to ob- 
tain spatially-resolved spectra of a scene, it can also oper- 



ate in an "occultation mode" where the short-wavelength 
VIS channel is turned off, while the longer-wavelength 
IR channel stares at a single pixel targeted at a star and 
obtains a series of rapidly-sampled near-infrared stellar 
spectra at 31 wavelengths between 0.85 and 5.0 microns. 
As the star moves behind the rings, its apparent bright- 
ness varies due to variations in the ring's opacity. Note 
that the response of the detector is highly linear, so af- 
ter a constant instrumental background is removed from 
each spectral channel, the data numbers returned by the 
instrument are proportional to the incident flux. In or- 
der to avoid contamination from sunlight scattered by 
the rings, we focus exclusively on data from one spec- 
tral channel covering the range 2.87-3.00 microns, where 
water ice is strongly absorbing. The rings are suffi- 
ciently dark at these wavelengths that ringshine is neg- 
ligible, and the measured signal is directly proportional 
to the transmission through the rings T. We can there- 
fore easily translate the raw data numbers into the slant 
optical depth along the line of sight through the rings 
T = — ln(T), or the normal optical depth of the rings 
T„ = Tsini? (i?=elevation angle of the star above the 
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rings). 

A precise time stamp is appended to each spectrum, 
facilitating the geometry reconstruction. Using the ap- 
propriate SPICE kernels, we can compute the radius 
where the starlight pierced the ringplane for each sam- 
ple in the given occultation. This calculation accounts 
for the light travel time from the rings to Cassini, and 
uses stellar positions taken from the Hipparcos catalog, 
corrected for parallax at Saturn. The positions of sharp 
edges of gaps and ringlets in the Cassini Division and C 
ring demonstrate that the resulting reconstructed geom- 
etry for each occultation is accurate to wit hin a kilometer 
(|Nicholson et al.ll2011l: [French et al.ll2011[ ). French et al. 
(in prep) have used a sub-set of these sharp features to 
make small corrections to the spacecraft's position dur- 
ing these occupations. These corrections not only reduce 
the dispersion in the radial position estimates of sharp 
edges to ^ 150 meters, but also yield more consistent 
estimates of the density wave phases (see below) . 

Between 2005 and 2009, VIMS obtained a total of 27 
occultation cuts through the C ring with sufhcient signal- 
to-noise and spatial resolution to discern the relevant 
waves. Table [5] provides a list of these occultations, iden- 
tifying the star observed, the "Rev number" (Cassini's 
orbit around Saturn) when the data were obtained, and 
whether the cut was obtained during ingress or egress. 
The table also provides the ephemeris time and observed 
inertial longitude when the star passed behind each of 
the six waves discussed below. Blank entries in this ta- 
ble correspond to cases where the occultation did not 
cover the particular wave or when a data gap corrupted 
the relevant profile, and thus are excluded from this anal- 
ysis. 

3.2. Initial examination of the waves using RS Cancri 
occultations 

Before describing the procedures we will use to deter- 
mine the m-numbers and pattern speeds of the waves 
from these occultation data, let us first examine some 
particularly informative occultations of the star RS Can- 
cri. Unlike the other occultations considered in this anal- 
ysis, these were chord occultations where the star cut 
obliquely through the rings. In such an occultation, the 
star passes through the same radial range at two very dif- 
ferent longitudes as the track enters and then leaves the 
rings. Comparisons of the wave profiles derived from the 
ingress and egress parts of these occultations thus pro- 
vide illustrative examples of how occultation data can 
constrain wave parameters. Furthermore, even simple 
visual inspections of these wave profiles indicates that 
several of the unidentified waves have very similar pat- 
tern speeds and likely the same m-numbers. 

Figures [2] and [3] show the relevant occultation profiles 
from the RS Cancri occultations observed during Revs 80 
and 85 (the Rev 87 RS Cancri occultation only probed 
waves exterior to 84,000 km and is not shown here). Of 
particular interest are the three waves W82.00, W82.06 
and W82.20, which are clustered within a region a cou- 
ple of hundred kilometers wide around 82,100 km. In the 
Rev 80 RS Cancri data the ingress and egress profiles for 
each of these waves are almost perfectly anti-correlated, 
with peaks in one profile corresponding to dips in the 
other and vice versa. This indicates that the phase differ- 
ence between these two occultation cuts is close to 180° 



for all three of these waves. This finding alone might be 
interpreted as a coincidence, but if we turn our attention 
to the Rev 85 data, we again find suspicious similarities 
among these three waves. For all three of the waves, the 
sharp peaks in the egress data occur about one-third of 
a cycle exterior to the sharp peaks in the ingress data. 
This again suggests that the phase difference between the 
ingress and egress cuts (in this case r^ 240°) is nearly the 
same for all three of these waves. 

This result is significant because it almost certainly 
implies that all three of these waves have the same m- 
number. The waves W82.00, W82.06, and W82.21 are 
close together in the rings. Hence, for each of these 
occultations, the ingress cuts for these three waves oc- 
curred at nearly the same longitude and time, and sim- 
ilarly the egress cuts are grouped closely in longitude 
and time (see Table [5]). This means the difference in the 
observed longitudes (SX) and the observed times (St) be- 
tween ingress and egress are almost identical for the three 
waves. Since the phase difference between the waves 
is given by S<j) — \m\{SX — flpSt), the easiest way for 
all three waves to have the same phase difference is for 
them to have the same 7n-number (which also implies, via 
Eq. [21 nearly identical pattern speeds) . Our full analysis 
of these waves confirms that this is indeed the case. 

Being able to determine that W82.00, W82.06 and 
W82.21 have the same tti- number from simple inspec- 
tion of the RS Cancri profiles is useful because this re- 
sult does not require the geometry of the occultations 
to be extremely accurate. The most likely errors in the 
occultation geometry correspond to shifts in the space- 
craft's position along its orbit, which would cause the 
ingress and egress occultation traces to shift in opposite 
directions in radius. Such errors will tend to shift the 
ingress profiles outward or inward relative to the egress 
curves by the same amount for all three waves. Although 
the ingress-egress phase differences would then change, as 
long as the wavelengths are similar this would not change 
the result that the phase differences between the ingress 
and egress cuts are nearly the same for all three waves. 

3.3. Computing phase differences between waves 

While simple inspection of the RS Cancri profiles sug- 
gests that several of the waves have similar pattern 
speeds and identical TTi-numbers, we need quantitative 
measurements of the phase differences between different 
occultation cuts to ascertain the actual values of m or Qp 
for these waves. Fortunately, we can estimate the phase 
difference between any two occultation cuts through a 
given wave using wavelet transforms. A wavelet is ba- 
sically a localized Fourier transform that can effectively 
cope with the rapid wavelength changes in typical spi- 
ral waves. Indeed, wavelet techniques have already been 
used to study density waves in order to obtain informa- 
tion about how the radial wavelength varies with distance 
from the resonance , which can constrain t he rings' sur- 
face inass density (JTiscareno et all l2007t iColwell et al.l 
I2009al : IBaillie et al.l 1201 II ). However, our analysis is dif- 
ferent in that it focuses on estimating the phase difference 
between two wave profiles rather than the wavelength or 
amplitude of the wave. 

We perform our wavelet analysis on the raw signal pro- 
files instead of the derived optical depth profiles shown in 
Figures [!]■ [3] In practice, this choice has little infiuence 
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Fig. 2. — Profiles of the relevant waves obtained during the occultation of the star RS Cancri on Rev 80. The black profiles were obtained 
during ingress, while the green profiles were obtained on egress. The normal optical depth values assume the star's elevation angle above 
the rings is 29.96°. Note that the phase differences between the waves seen in ingress and egress are the same for the three waves found 
around 82,000 km. This suggests that all these waves have the same m-numbers. 



on the derived phase differences (the conversion from T 
to T just introduces an overall constant phase shift of 
180°), but it does eliminate any possible complications 
that could arise due to the uncertainties in the unoc- 
culted star signal. Before applying the wavelet transform 
to the relevant occultation data, we first interpolate the 
data onto a uniform array of radii with a radial spacing 
of 50 m. A small part of this re-sampled profile centered 
on each wave (see Table |3] for the exact ranges) is then 
fed into the publicly-available IDL routine wavelet (see 
Torrence and Compo 1998, the default Morlet mother 
wavelet with wq = 6 is used throughout this analysis). 
The resulting wavelet is a two-dimensional array of com- 
plex numbers as a function of radius and radial wavenum- 
ber VV(r, /c) . (Note we use the wavenumber k here instead 
of the wavelength A in order to avoid any possible confu- 
sion with the radiation's wavelength or the longitudinal 
coordinate.) Let us denote the real and imaginary parts 
of the wavelet as Wr and W/ , respectively. We can then 
define the wavelet power as Vir, k) = W^ + VVf and the 
wavelet phase as (p{r, k) = tan~^(W7, Wr). 

Figure m illustrates how the wavelet phase and power 
vary with position and wavelength across the wave. As 



expected there is a diagonal ridge in the wavelet power 
that tracks the observed trends in the wave's wavelength. 
(Recall that for a wave generated at an OLR, the radial 
wavelength decreases as the wave propagates inwards.) 
Furthermore, we can observe that where the wave is 
strong and the wavelet power is high, the contours of 
constant phase are nearly vertical, so a well-defined phase 
can be ascribed to every radius in the wave. The values 
of the phase in this region are also reasonable given the 
profile, being near ±180° at local minima and around at 
local maximao Thus, for such waves we can reduce the 
two-dimensional wavelet to estimates of the wave power 
and phase as a function of radius by appropriately aver- 
aging over a range of spatial wavenumbers. 

In order to filter out large scale background variations 
in the rings' opacity while still capturing most of the 
wave's power, we include a limited range of wavenumbers 
in these averages. Specifically we exclude all wavenum- 
bers less than fci = 27r/(5 km) and all wavenumbers 
greater than fc2 — 27r/(0.1 km). The effective power of 

■^ The wavelet phase also increases with increasing radius, con- 
sistent with the expected trend in (/'(r) for a density wave. 
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Fig. 3. — Profiles of the relevant waves obtained during the occultation of the star RS Cancri on Rev 85. The black profiles were obtained 
during ingress, while the green profiles were obtained on egress. Note that the phase differences between the waves seen in ingress and 
egress are again the same for the three waves found around 82,000 km. This implies that all these waves have the same m-numbers. The 
two waves found outside 84,000 km also show similar phase differences, indicating that they may have same m-number as each other. 



the wave at a given radius is therefore defined to be: 



P,s{r)^MY.r{r,k), 



(8) 



k=ki 



where AT is a normalization constant. Note that since 
we only are interested in relative power for this analysis, 
we choose TV so that PcS is equal to unity at its peak 
value. Also note that this sum is done over a series of 
logarithmically-spaced values of k. 

In order to derive an effective phase for the wave at 
each radius, we first compute an effective average real 
part and imaginary part of the wavelet: 



WRj{r) 



j:^B.j{r,k)V{r,k) 
Y.V{r,k) 



(9) 



where the sums are again over all values of k between fci 
and fc2- Note these averages are weighted by the power, 
so that the regions with the strongest signals dominate 
the averages. From these average wavelet components, 
we can compute the average phase at each radius: 



i&n-\Wi,WR) 



(10) 



Note by computing the average of components W/ and 
Wr instead of averaging the local phases (f{r,k), we 
avoid any difficulties involved in averaging a cyclic quan- 
tity. The phase difference between two occultations 
should then simply be the difference in the two values 
of (/'(r). 

Figures [5] and [5] illustrate the results of these calcu- 
lations for two of the waves seen in the RS Cancri oc- 
cultation from Rev 85. As expected, there is a peak in 
the wave power Pcff located near the center of the wave. 
Also, the difference in the wave phase between the egress 
and ingress cuts is relatively constant in the region where 
the wave signal is evident, as desired. Furthermore, the 
numerical values of these phase differences are consistent 
with the observed profiles. For the W84.64 wave in Fig- 
ure [6l the phase difference is around 180°, which is what 
one would expect given that the peaks in the ingress pro- 
file occur at the same location as the dips in the egress 
profile and vice versa. Similarly, the W82.21 data in Fig- 
ure [5] yield a phase difference of around 240° , which is 
consistent with the dips in the egress profile always lying 
about 1/3 of a cycle exterior to those in the ingress pro- 
file. This gives us some confidence that our rather simple 
approach can extract useful phase information from these 
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Fig. 4. — Plot showing the wavelet amplitude and phase derived from the W84.64 wave observed by the Rev 106 RCas observation. The 
top panel shows the occultation profile (in raw Data Numbers, which is proportional to transmission) as a function of radius. The bottom 
panel shows the wavelet phase and power as functions of radius and spatial wavelength. The wavelet phase is indicated by greyscale levels 
(black= —180°, white= +180°) while the overlaid green contours are levels of constant wavelet power. The peak wavelet power follows a 
diagonal ridge that corresponds to the wave's increasing wavelength with radius. Note that where the wavelet power is strong, the contours 
of wavelet phase are nearly vertical and correspond to the expected phase of the wave {e.g., the phase wraps from -180° to 180° at locations 
corresponding to sharp minima in the profile). 



wave profiles. 

For the purposes of this analysis, we require a single 
estimate of the average phase difference for any pair of 
occultation cuts. We estimate this parameter as the 
weighted average of the phase differences S(f>{r), with 
each of the individual values weighted by the average of 
the two Peff curves in order to ensure that regions with 
high signal contribute most to the final estimate. In or- 
der to further reduce the possibility of contamination in 
the phase difference estimate, we only consider regions 
where the average power of the two waves is more than 
0.9 (this threshold is indicated by a dotted line in Fig- 
ures [5] and [ni and the averaging region is demarcated by 
vertical dashed lines). Due to the weighting, the result- 
ing estimates of the phase differences are not particularly 
sensitive to the exact value of this threshold. In addition 
to computing the weighted average phase difference, Scj) 
between the two cuts (which is shown as a dotted line in 
the bottom panel of Figures [5] and [5]) , we also compute 
the standard deviation of the phase difference values in 
the selected region at/,, which quantifies the reliability of 
the phase estimate. Note, however, that a^ will under- 
estimate the uncertainty in the 6(j) because it does not 
include the efii'ects of uncertainties in the reconstructed 
geometry. Such geometric uncertainties are difficult to 
quantify a priori^ and will be considered in more detail 
below. 

To recap this procedure, the steps are: (1) Compute 
the wavelet transform, W{r, k) for each observation and 
the corresponding power spectrum, V{r,k). (2) Com- 



pute weighted averages for Wii{r) and Wi{r), and thence 
the average phase 0(r). (3) Compute the average radial 
power profile, Pcs{r). (4) Compute the average phase 
difference between two wave profiles (50, weighted by the 
average of the two power profiles. 

3.4. Using phase differences to constrain m-numbers 
and pattern speeds 

In order to illustrate how quantitative estimates of Scj) 
can be used to constrain ?7T,-numbers and pattern speeds, 
let us examine the phase differences between the various 
RS Cancri ingress and egress profiles listed in Table [2j 
Note that the phase differences for the W82.00, W82.06 
and W82.21 waves are indeed very similar for both the 
Rev 80 cuts (being around ~ 150°), and the Rev 85 cuts 
(where they are all ~ 250°). This is consistent with the 
above visual inspection of the profiles and supports our 
contention that these three waves must share a common 
m-number. However, with these measured values of 6(j), 
we can now also identify which values of m and flp come 
closest to predicting the observed phase differences. 

For any given wave, the local mean motion n and ap- 
sidal precession rate w of the ring material are straight- 
forward functions of the wave's resonant radius r^ and 
Saturn's gravitational field ( Murray fc D ermott 1999). 
Hence, i f we assume th e grav itational field parameters 
given in iJacobson et al.l ()2006D . we can use Equation [5] 
to compute the expected pattern speed l^p for any given 
value of m at any of the relevant resonant radii. Fur- 
thermore, for a particular pair of occultation cuts, the 
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TABLE 2 

Phase differences between ingress and egress cuts for the RS Cancri occultations 



Wave 


Resonant 


Region 


Rev 80 


Rev 85 


Rev 87 


Possible 




Location" 


Considered'' 


50 


5<f> 


&(t> 


m-valucs'= 


W80.98 


80988 km 


80970-80995 km 


141.6° 


— 


— 


-4, -3, +6 


W82.00 


82010 km 


81992-82012 km 


148.5° 


250.6° 


— 


-3, +6 


W82.06 


82061 km 


82040-82065 km 


148.1° 


249.1° 


— 


-3, +6 


W82.21 


82209 km 


82190-82215 km 


150.1° 


251.1° 


— 


-2, -3, -f 6 


W84.64 


84644 km 


84625-84650 km 


145.4° 


195.5° 


235.5° 


-2, 4-5, -1-6 


W87.19 


87189 km 


87175-87205 km 


d 


(182.2°)"= 


191.4° 


-2, -H5 



" Inferred resonance locations from iBaillie et aP l|20111 ) Table 7 

* Region over which wavelet power and phase are computed. 

° m-values between -10 and +10 where the predicted &4> values are within ±30° of the observed values. 

"* Phase difference not determined due to data gap 

*= Phase difference computed using wavcvectors between k = 27r/(0.1 km) and k = 27r/(2 km) because otherwise combined peak average 

power is never above threshold. 



Ingress 



Rev 85 RSCnc 




Ingress 



Rev 85 RSCnc 



Egress 



82.195 82.200 82.205 82.210 
Radius (1000 km) 

Fig. 5. — Results of the wavelet calculations of the phase differ- 
ence in wave W82.21 between the ingress and egress cuts from the 
Rev 85 RSCnc occultation. The top panel shows the two occulta- 
tion profiles, while the middle panel shows the integrated wave 
power Pcff('') between wavenumbers of 27r/(5 km) and 27r/(0.1 
km). The bottom panel shows the phase difference Srf>{r) = 
4'egress{'>') — 4'ingreas{r) between these two cuts (see text for ex- 
planations of the dashed and dotted lines) . Note that the average 
phase difference is computed using only the data where the average 
Pelf of the two profiles is above 0.9. The average phase difference 
is near 240° , which is consistent with the offset between the ingress 
and egress wave profiles noted in Fig.[3l 



difTerence in the observed ring longitudes SX and the dif- 
ference in observation times St can be obtained from Ta- 
ble m We can therefore use the computed pattern speed 
to calculate the expected phase difference between any 
two occultation cuts (S0(predicted) = |?ti|(JA — fipSt) for 
each value of m (cf. Eq. [7]). 

Figures [7] and [S] show the differences between the ob- 
served and predicted values of 6(j) for the waves W82.21 




84.625 84.630 84.635 84.640 84.645 84.650 
Radius (1000 km) 

Fig. 6. — Results of the wavelet calculations of the phase dif- 
ference in wave W84.64 between the ingress and egress cuts from 
the Rev 85 RSCancri occultation, following the same layout as in 
Fig-El I^ this case, the average phase difference is near 180°, which 
is again consistent with the ingress and egress wave profiles. 

and W84.64 as a function of mQ For any given observa- 
tion, the difference (5(/)(observed) - 6(j) (predicted) cycles 
repeatedly through 360° with increasing or decreasing m. 
Hence there are multiple possible m-values that could 
be consistent with the single observed phase difference. 
However, different observations have different values of 
St and 6X and thus show different trends in these plots. 
Hence there are relatively few m-values that could be 
consistent with the two or three observations illustrated 
in Figures [7] and m If we conservatively assume a ±30° 

•^ Note that the time elapsed between the ingress and egress cuts 
can be as much as two hours. This is a non- trivial fraction of the 
6.5-7 hour orbital period in this region and therefore cannot be 
ignored in the calculation of the predicted S<p. 
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Fig. 7. — Plot showing the difference between the observed and predicted values of Sep between the ingress and egress cuts of the RS 
Cancri occupations for the wave W82.21 as a function of the assumed m-number, given the stipulated S\ and St values. Different symbols 
correspond to different pairs of occultation cuts. Note that the residuals for both observations are close to zero when m = —2, —3 and +6, 
so these values of m are the ones most consistent with the observed phase differences. 



W84.64, Assumed r^ = 84644.0 km 

Rev 80 RSCnc, d,p = 1 45.4°, (5t = 71 71 .9 sec, dX = 45.8°, 

Rev 85 RSCnc, ^0 = 1 95.5°, (5t = 5298.9 sec, (5A = 33,3°, 

Rev 87 RSCnc, d,p = 235.5°, dt = 4146.5 sec, dX = 25.9° 




Fig. 8. — Plot showing the difference between the observed and predicted values of 50 between the ingress and egress cuts of the RS 
Cancri occupations for the wave W84.64 as a function of the assumed m-number. Note that the residuals for all three observations are 
close to zero when m = —2, +5 and +6, so these values of m are the ones most consistent with the observed phase differences. 
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uncertainty in the phase determinations (which is con- 
sistent with the results of the global analysis described 
below), then to = -3 -2 and +6 are consistent with the 
RS Cancri observations of the W82.21 wave, and m = - 
2, +5 and +6 are consistent with the observations of the 
W84.64. The other waves yield similar results, as shown 
in Tabled! Note that the W82.00, W82.06 and W82.21 
waves are consistent with the same set of TO,-values, as 
expected 

The small number of to,- values that are consistent with 
this limited number of 6(f> estimates demonstrates that 
our measurements of Scj) can constrain the symmetry 
properties and pattern speeds of these waves. However, 
it is also clear that we cannot uniquely determine the 
771- value for any of these waves with only the RS Cancri 
data. Thus the next step in this analysis is to extend this 
approach to include all pairs of occultation cuts listed in 
Table p 

While the procedures described above allow us to cal- 
culate 5(j) for any possible pair of occultation cuts listed 
in Table m in practice some occultation pairs do not yield 
reliable phase difference estimates. For some pairs, the 
average wavelet power for the two cuts never exceeds 0.9 
of the peak power in each cut. This implies that the 
peak wave signal in the two cuts is coming from different 
parts of the wave (perhaps because they have different 
intrinsic resolutions), and so we do not consider such 
pairs. We also exclude any occultation pairs which yield 
a rms phase-difference scatter a^ greater than 20° be- 
cause large values of a^ indicate that the wavelet was 
unable to identify a consistent phase difference (perhaps 
due to small gaps or cosmic rays in one of the profiles) . 
Finally, we only consider pairs of occultations where the 
time difference is less than 300 days. This prevents any 
aliasing that might occur due to the limited number of 
observations with larger time separations. After apply- 
ing these selection criteria, we have between 100 and 260 
6(1) estimates for each of the waves (see Table [3] and Ta- 
ble E]). 

As with the RS Cancri data described above, we con- 
strain the symmetry properties and pattern speeds of 
each wave by computing expected values of 6(1) for dif- 
ferent combinations of m and Up and comparing these 
numbers with the observed 6(j) values. Specifically, we 
seek values of m and Vlp that minimize the rms resdiuals 
of 5(/)(observed) - 6(j) (predicted). In principle, we could 
just compute these rms variations for each to, assuming 
ftp is given by Equation [21 as we did for the RS Cancri 
data. However, given the uncertainties in the planet's 
gravitational field and the precise locations of the reso- 
nant radii for these waves, we will instead compute the 
rms phase difference residuals over a finite range of pat- 
tern speeds surrounding the expected ftp corresponding 
to each m. 

Since we are studying patterns that rotate around the 
planet at hundreds of degrees per day, and we are us- 
ing occultations separated in time by up to 300 days, it 
is important to verify that our 6(j) calculations are suffi- 
ciently accurate to yield meaningful constraints on m and 
Qp. This is most easily done by first testing our methods 
on waves with known pattern speeds. We therefore ap- 
plied the above algorithms to occultation profiles of the 
Prometheus 8:7 density wave in the outer A ring. The top 
panel in Figure [9] shows the rms phase difference resid- 



ual (observed-predicted) as a function of pattern speed, 
assuming the pattern has the expected m = 8. Note the 
sharp dip at 587.30°/day, which corresponds almost ex- 
actly to Prometheus' mean motion and the expected pat- 
tern speed for this wave. Thus the pattern speed with the 
minimum variance in the phase residuals matches the ex- 
pected perturbation frequency, as desired. Furthermore, 
no similar dip is seen when other values of to are tried. 
These results demonstrate that the above procedures can 
indeed yield reliable estimates for the pattern speeds and 
?Ti-numbers of even tightly-wound density waves. 

The lower panel of Figure |9] shows the phase differ- 
ence residuals as a function of time separation 6t. These 
points are randomly scattered about zero, as desired, but 
with a ±20° scatter which probably reflects not only sta- 
tistical uncertainties in the phase estimates themselves 
but also uncorrected systematic errors in the occultation 
geometries. As noted above, the latter can shift wave 
profiles relative to each other, which will produce errors 
in the phase differences that are hard to model. However, 
these data demonstrate that these systematic errors do 
not prevent us from determining the pattern speeds of 
spiral waves with wavelengths that are sufficiently long. 
More specifically, the above procedures should yield reli- 
able results so long as the uncertainties in the occultation 
geometry do not produce phase errors exceeding 90° . As 
discussed in Section 12.31 above, this should be the case 
for the six waves considered here. 

The small dispersion in the phase difference residu- 
als not only requires accurate phase difference measure- 
ments, the phase differences must also follow the pre- 
dicted trend given by Eq. [7] over a sufficiently long pe- 
riod of time. For waves generated by satellites, the latter 
condition is equivalent to assuming the relevant moon 
has a constant mean motion, which is reasonable for 
most of Saturn's moons (the obvious exceptions being 
the co-orbitals Janus and Epimetheus, see Tiscareno et 
al. 2006). For waves generated by planetary normal 
modes, this condition requires assuming that the relevant 
oscillation has a constant frequency and a coherent phase 
over the time period spanned by the observations. This 
assumption not only turns out to be consistent with the 
observations (see below), it can also be justified a priori 
based on considerations of the waves' group velocity. 

Density waves in dense rings propagate away from the 
resonance at a finite speed, so any shift in the frequency 
or phase of the perturbing potential generates dis continu- 
ities in the wave's profile (jTiscareno et al.ll2006( ). These 
features propagate through the rings at the group veloc- 
ity Vg ~ nGao/K, where G is the universal gravitational 
constant, ctq the ring's average surface mass density, and 
K the radial epicyclic frequency (Shu 1984). For the rele- 
vant C-ring waves, ctq is between 1 and 10 g/cm'^ (Zebker 
et al. 1985, Hedman et al. 2011, Section [S] below) , and 
Vg is between 0.4 and 2.5 km/year. The relevant waves 
extend 10-20 km from the resonances and they do not ex- 
hibit any deviations from the expected smooth trends in 
their wavelengths. Thus, regardless of their origins, the 
frequencies and phases of the relevant perturbing poten- 
tials appear to be coherent over several years. Hence we 
are justified in assuming that the phase differences can 
be predicted using Eq. [7] for time intervals of one year or 
less. 
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Prometheus 8:7 Wave, Assumed Resonant Radius = 127613.5 km, Assumed m = 8 
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Fig. 9. — A test of our pattern-speed determination algorithms using the Prometheus 8:7 wave in the inner A ring. The top panel shows 
the rms phase difference residuals as a function of pattern speed assuming the pattern has an m = 8, as expected for this wave. The dashed 
line marks the predicted pattern speed for this pattern, while the dotted line marks the pattern speed that gives the minimum variance (in 
this case, these two lines are almost on top of each other). The bottom panel shows the residuals in the phase differences from this best- fit 
solution as a function of time difference between the pairs of observations. The scatter in these data likely represents residual geometrical 
uncertainties in the various profiles. 



4. RESULTS 

We searched for the best-fitting combinations of m- 
numbers and pattern speeds for all six of the waves 
described in Section 12.31 Given that the wavelength 
of these waves increases with increasing radius, we ex- 
pected these patterns would be generated by outer Lind- 
blad resonances, and so the 77i-numbers of these waves 
wo uld be negat i ve. A lso, previous studies of t hese waves 
bv lRosen et all (|1991D and lBailhe et al.l (|2011l ) suggested 
that if |?Ti| was above 5, the implied surface mass density 
of the C ring would be unreasonably high (see Section [5] 
below). Finally, the comparison of ingress and egress 
phases for the RS Cancri occultations in Table [2] suggests 
plausible values for m of —2,-3 and —4. Thus we ex- 
pected that for these waves m would lie between —2 and 
—5. Nevertheless, for the sake of completeness, we con- 
sidered all m values between -1-10 and —10, and searched 
for minima within 10° /day of the expected pattern speed 
given by Equation [50 

For each wave, we found a strong dip in the residual 
variance for only a single value of m. For wave W80.98, 
this dip occurred with m = —4; for each of waves W82.00, 
W82.06 and W82.21 it happened with m = -3; and 
for waves W84.64 and W87.19 it was m = -2. Fig- 
ure [TU] shows profiles of the rms residuals versus pat- 
tern speed for all six waves, assuming the appropriate 

* We also looked for patterns that would be produced by the 
higher-order Enceladus 5:1 (m = 2, Qp = 656.7° /day) and the 
Janus 5:2 (m = 3, Qp = 863.7° /day) inner vertical resonances 
that lay close to the W80.98 and W82.06 waves. No significant 
dips in rms residuals were found. 



771- value. The resulting best-fitting pattern speeds can be 
compared with those calculated for each wa ve, using the 
resonance locations specified in Table 7 of iBaillie et al.l 
(|201lD . For all of the waves the minimum in residuals oc- 
curs very close to the expected pattern speed. For five of 
the waves, the best-fit pattern s peed is within 0.08° /day 
of the value predicted using the IBaillie et al.l (|2011f ) res- 
onance locations, while for the W87.19 wave the best- fit 
pattern s peed is 0.16°/day lo wer than one would predict 
using the IBaillie et al.l ()2011h resonance location. These 
small offsets could represent small errors in either the as- 
sumed gravity field or the estimated resonant radii. For 
example, the offset between the best-fit and predicted 
pattern speeds for the W87.19 wave could be resolved if 
the real reso n ant ra dius r^ lies just 4 km exterior to the 
IBaillie et ahl ()2011[ ) estimate of 87,189 km. This wave 
is superimposed on a peak in optical depth (minimum 
in transparency), and the resulting variations in the sur- 
face mass density across the wave may have complicated 
earlier efforts to estimate the resonant radius. Indeed, 
looking at Figures [U [2] and [31 the outer edge of this 
wave could easily fall somewhere around 87,193 km. 

Figure [TT] shows the individual phase difference resid- 
uals as a function of St for the best-fit solutions. These 
residuals typically have a spread of approximately ±45°, 
which is noticeably larger than that found for the 
Prometheus 8:7 wave discussed above. This is not en- 
tirely surprising given the generally shorter wavelengths 
and more limited extents of the C-ring waves. Fewer 
cycles are visible in most occultation profiles of the C- 
ring waves than for the A-ring wave, so the estimates of 
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Fig. 10. — Plots showing the rms phase difference residuals as a function of pattern speed for each of the six waves, assuming the pattern 
has the indicated m -numljers. The dashed line marks the predicted pattern speed for this pattern at the resonant location provided by 
IBaillie et al.l 1)20111) . while the dotted line is the pattern speed that gives the minimum variance. 



the phase differences should be more uncertain. Further- 
more, the maximum wavelength of the Prometheus 8:7 
wave exceeds 10 kilometers, while the C-ring waves all 
have maximum wavelengths that are less than 5 kilome- 
ters. The estimated phases of the C-ring waves should 
therefore be more sensitive to small errors in the geo- 
metrical reconstructions. Indeed, if we do not include 
the small corrections to the occultation geometry pro- 
vided by French et al. (in prep.), the rms residuals of 
the best fit solutions increase by roughly a factor of two. 
However, even with this increased dispersion, the best- 
fit solution for each of the six waves remains close to the 
predicted pattern speed for the appropriate OLR. Hence, 



our identification of the best-fit m-numbers and pattern 
speeds appears to be robust against any residual uncer- 
tainties in the reconstruction geometry. 

5. DISCUSSION 

Table [3] summarizes the results of this analysis, includ- 
ing the TO- numbers, pattern speeds ftp, and pattern rota- 
tion periods {2n/flp) for the six waves. We may conser- 
vatively estimate the uncertainties in the pattern speeds 
to be 0.5°/day, since an error in the pattern speed of this 
magnitude would produce a trend in the residuals that 
would be clearly detectable in the data (the residuals 
would reach ~180° at 300 days). This is also compa- 
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Fig. 11. — Plots showing phase difference residuals (observed-predicted) for each of the six wave s, assuming each pattern has the indicated 
m-numbers and the pattern speeds, which correspond to the best-fit values shown in Figure [TOl 

rable to the widths of the minima in Figure [TU] While 
the relatively tight distribution of the residuals already 
gives us some confidence that we have determined the 
correct values of m and fip for each of these waves, the 
following considerations lend additional support to these 
identifications: 

1. Each of the pattern speeds is consistent with that 
calculated for an outer Lindblad resonance at that 
location, for the appropriate m-number. 

2. Our full analysis indicates that the waves W82.00, 
W82.06 and W82.21 all have m = -3. This is 
consistent with the similar observed phase shifts 
between ingress and egress observed for these three 
waves in the RS Cancri occulations discussed in 



Sections 13.21 and 13.41 Recall that these arguments 
were robust against uncertainties in the occultation 
geometry, and thus provide a useful check on the 
above analysis. 

3. The derived values of m yield plausible values for 
the C-ring's surface mass density. 



The last of these 
tion and qualificat 
wave's wavelength 
of the ring's averaf 
its mass extinction 
ever, the estimated 
m, so prior to this 
estimates oi ao/\m 



three points requires some explana- 
lon. In general, the variation of a 
with radius can provide an estimate 
z,e local surface mass density ao and 

coefficient Tu/cto, via Eq. ^. How- 
values of these parameters depend on 

work researchers could only provide 
— ll for these six unidentified waves. 
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TABLE 3 

Results of the present analysis of the six waves 



Wave 


Resonant 


Region 


Af((50)= 


m 


Pattern 


Rotation 




Location" 


Considered'' 






Speed"^ 


Period 


W80.98 


80988 km 


80970-80995 km 


136 


-4 


1660.3° /day 


312.2 min 


W82.00 


82010 km 


81992-82012 km 


219 


-3 


1736.6° /day 


298.5 min 


W82.06 


82061 km 


82040-82065 km 


217 


-3 


1735.0° /day 


298.8 min 


W82.21 


82209 km 


82190-82215 km 


191 


-3 


1730.3° /day 


299.6 min 


W84.64 


84644 km 


84625-84650 km 


257 


-2 


1860.8° /day 


278.6 min 


W87.19 


87189"= km 


87175-87205 km 


111 


_2 


1779.5°/day 


291.3 min 




ible 7 










n|Baillieetal.H2011|)T 





* Region over which the wavelet phase and power are computed. 

'^ Number of Sep estimates used in the fits. 

"* Best-fit value. A conservative estimate of the uncertainties in these numbers is 0.5° /day, which would correspond to a ~ 180° residual 

at 300 days. 

° The true value may be closer to 87193 km (see text). 



Table m list s the estimates oi ao /\m — 1| for all six waves 
derived by iBaillie et al.l ()2011[ ). and the resulting esti- 
mates of (Jo and t„/(To assuming the m-numbers derived 
here0 These values can be compared with estimates de- 
rived from other features in the rings, but in making 
these comparisons we must keep in mind that both cto 
and T„/(To may vary with position across the ring. 

Studies o f the five identifiable C-ring waves by 
IBaillie et al.l ()2011[ ) yield extinction coefficients ranging 
between 0.13 and 0.36 cm^/g, which are higher than our 
estimates. However, most of these identifiable waves oc- 
cur in very different environments from those occupied 
by the six unidentified waves considered here. The two 
waves found near the Mimas 6:2 and Prometheus 4:2 
resonances occur within a high-opacity plateau, while 
the waves associated with the Titan nodal resonance and 
the Mimas 4:1 resonance are found in rather low optical 
depth regions interior to 78,000 km. This leaves only the 
wave associated with the Atlas 2:1 ILR, which is within 
500 km of W87.19 and not in a re gion of obviously el - 
evated optical depth. Intriguinglv. IBaillie et al.l ()2011l ) 
estimate that t„/(7o ~ 0.19 ± 0.04 cm^/g for the Atlas 
2:1 wave, which is not too different from the value of 0.11 
cm^/g derived above for the nearby W87.19 wave. 

For the waves between 79,000 km and 85,000 km there 
are no known density waves that can provide independent 
estimates for ctq or T„/cro, but estimates of these param- 
eters have been extracted from other types of observa- 
tions. Based on the particle size distribu tions derived 
from V oyager RSS radio occultation data, iZebker et al.l 
(|1985| ) estimated the average surface mass density of the 
C ring between 78,429 km and 84,462 km to be 3.2 ± 1.8 
g/cm^. More recently, iHedman et al.l ()2011|) found that 
variations in the wavelength of a vertical corrugation ex- 
tending between 78,000 and 84,000 km indicated that 
this part of the ring has a mass extinction coefficient of 
roughly 0.02 cm^/g. Both of these numbers are reason- 
ably consistent with the values of ctq and t„/(To given in 
Table 0] for the inner five waves. We thus conclude that 
the 771- values derived above do indeed yield plausible val- 
ues for the mass density of the middle C ring. 

Assuming the above rn-numbers and pattern speeds 
are correct, we may now ask what could be producing 
these waves. They almost certainly cannot be produced 
by any of Saturn's moons, since such a moon would need 

^ Note that, for an OLR, |m — 1[ becomes \'m\ + 1. 



to be orbiting just above Saturn's cloud tops, or within 
the tenuous D ring. Not only have no such moons been 
detected by Cassini, despite extensive imaging of this 
region, but any km-size icy or rocky objects unfortunate 
enough to find themselves here would be quickly torn 
apart by tidal forces. 

Instead, these waves may be generated by perturba- 
tions from oscillating normal modes within the planet. 
Indeed, using a particular model f or Saturn's inter- 
nal structure, IMarlev fc Pored (|1993[ ) predicted pattern 
speeds for a number of normal modes that could pro- 
duce several of the w aves identified in the Voyager radio 
occultation (jRosen e t al. 1991,). Those authors even sug- 
gested that W80.98 might be an w = —4 wave generated 
by a resonance with the I = 4 sectoral {i.e., m = I) f- 
mode in the planet (which they denote 4/), and that 
either W82.06 or W82.20 might be generated by a sim- 
ilar I — m ^ S planetary /-mode. Indeed, our fitted 
pattern speeds appear to be very consistent with Marley 
and Porco's predictions for the pattern speeds generated 
by normal mode oscillations inside the planet. Their pre- 
dicted rotation period for the 4/ mode of 310.6 minutes 
is quite close to the observed period of 312.2 minutes for 
the m — —4 wave W80.98. Their predicted 3/ mode 
period of 296.4 minutes is close to the periods of the 
m = -3 waves W82.00, W82.06 and W82.21, which are 
298.5 — 299.6 minutes. Finally, the predicted 2/ mode 
period of 286.4 minutes lies in between the observed pe- 
riods of our two m = —2 waves W84.64 and W87.19. 
Hence, it seems quite likely that /-mode acoustic oscilla- 
tions within the planet may be producing these six waves, 
in which case the precise pattern speeds presented here 
should be able to help constrain interior models for Sat- 
urn. (As but one example, we note that , with in the range 
of pattern speeds calculated bv lMarlevI ()1991D for various 
Saturn models, those which best fit the observed wave lo- 
cations correspond to models with no interior differential 
rotation.) 

However, while each of the six waves has a pattern 
sp eed that fits with one of the normal modes studied 
by IMarlev fc Porcol (|1993f) , the number of waves associ- 
ated wi th each normal mode i s much harder to under- 
stand. IMarlev fc Porcol (|1993D predicted that a single 
wave would be associated with each normal mode oscil- 
lation. Furthermore, only the /-modes with m = I were 
predicted to have resonances in the C-ring (modes with 
/ > TO or with internal radial nodes were predicted to 



16 



TABLE 4 
Ring mass density and extinction coefficients 



Wave 


m 


<To/|m-l|" 


r„" 


(T„/ao)|m-ir 


o-() 


r„/o- 


W80.98 


-4 


1.17 g/cm^ 


0.13 


0.11 cm^/g 


5.85 g/cm^ 


0.022 cm^/g 


W82.00 


-3 


1.42 g/cm^ 


0.14 


0.10 cm^/g 


5.68 g/cm2 


0.025 cm2/g 


W82.06 


-3 


2.54 g/cm2 


0.28 


0.11 cm^/g 


10.16 g/cm2 


0.028 cm2/g 


W82.21 


-3 


1.73 g/cm2 


0.13 


0.08 cm2/g 


6.92 g/cm2 


0.020 cm2/g 


W84.64 


-2 


1.35 g/cm2 


0.11 


0.08 cm^/g 


4.05 g/cm^ 


0.027 cm^/g 


W87.19 


-2 


0.47 g/cm2 


0.15 


0.33 cm2/g 


1.41 g/cm^ 


0.11 cm^/g 



" from Table 7 of lBaillie et"an ||20T!1 '). 

have significantly faster pattern speeds and mucli smaller 
gravitational signatures). Thus they expected to find a 
single wave with m — —2, —3, —4, etc. Instead, we find 
three m — —3 waves in close proximity to each other, 
and two m = —2 waves with a larger separation. This 
suggests that Saturn's internal structure is more complex 
than Marley and Porco assumed. Thus future efforts to 
interpret these findings will likely need to consider the 
effects of differential rotation, compositional gradients, 
or a solid core. 

At the same time, further investigations of the C-ring 
waves should provide additional constraints on poten- 
tial models of Saturn's interior. For example, once more 
accurate geometric reconstructions of the stellar occul- 
tations are available, it should be possible to examine 
the pattern speeds of shorter- wavele ngth waves, such as 
the waves designated "d" and "h" bvlRosen et al] ()1991f ) 
or W85.67 and W83.63 bv lColweh et al.l (I2009bn . These 
particular features coul d represent add itional OLRs with 
|to| > 4 predicted bv IMarlev fc Porco (1993) . Mean- 
while, structures interior to 79,000 km could be gener- 
ated by vertical resonances with a different class of plan- 
etary normal modes. There could even be additional 
m = —4, —3 or —2 patterns lurking within the rings that 
are either too subtle to detect in individual profiles or 
obscured by other ring features. 

These waves also potentially contain information about 
the amplitude of the planetary oscillations. All else being 
equal, a larger-amplitude oscillation in the planet should 
generate larger fractional density variations in the cor- 
responding density wave. This would imply that the 
3/ osciUations responsible for the W82.21 and W82.06 
waves are larger than the one generating the W82.00 
wave, and that the 2/ oscillation responsible for W84.64 
is larger than the one that gives rise to W87.19. How- 
ever, the amplitude of the wave depends not only on 
the driving torque, but also on how quickly collisions 
among the ring particles dis sipate coherent motions (jShul 
119841: iTiscarenoet aLll2007[ ). This comphcates any ef- 
fort to extract quantitative estimates of the perturbing 
force from the wave amplitude, especially for waves that 
exhibit large fractional density variations, like the ones 
considered here. 



Finally, the ring data can constrain how quickly the 
oscillations are generated and dissipated within Saturn's 
interior. As mentioned in Section [3.41 above, the smooth 
trends in the waves' wavelength with radius indicate 
that each perturbing potential has maintained a coherent 
phase for several years. The small dispersion in the phase 
residuals illustrated in Figure [TT] confirms this supposi- 
tion, and furthermore we find no statistically significant 
trends in the dispersion with St out beeyond 300 days. 
If we assume that the planetary oscillations are stochas- 
tically excited, then this long coherence time implies a 
correspondingly high quality factor Q for these oscilla- 
tion modes. Indeed, given that the oscillation periods 
of the re levant modes in the p lanet's frame are ~ 200 
minutes (jMarlev fc Porcolll993fl . a coherence time in ex- 
cess of 300 days would imply a Q > 10,000. This num- 
ber is comparable to limits computed from the expected 
tidal evolut ion of Saturn's moon s over the age of the so- 
lar system (jDermott et al.l 11988^. and is comparable to 
estimates of Jupiter's Q deri ved from astrometric satel- 
lite data (jLainev et al.ll2009l ). However, recent analyses 
of astrometric measurements of Saturn's m oons indicate 
that Saturn's Q is between 1000 and 2000 (jLainev et al.l 
l2012f ). It is not obvious how this surprisingly low es- 
timate can be reconciled with the long coherence times 
observed here. Nevertheless, further analyses of the full 
span of Cassini occultation data should provide novel 
constraints on dissipation rates within Saturn. 
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TABLE 5 
Observed times (in seconds of ephemeris time, measured from the J2000 epoch) and inertial longitudes (measured relative 

TO the longitude of ascending node on J2000) FOR THE VARIOUS OCCULTATION CUTS THROUGH EACH WAVE. 



Star 


Rev 


i/c 


W80.98 


W82.00 


W82.06 


W82.21 


W84.64 


W87.19 


RHya 


036 




220948260. 
174.072° 


220948093. 
175.060° 


220948085. 
175.108° 


220948061. 

175.249° 


220947679. 
177.407° 


220947295. 
179.453° 


a Aur 


041 




227949007. 
343.813° 


227948789. 
345.077° 


227948779. 
345.137° 


227948748. 
345.316° 


227948261. 
348.013° 


227947783. 
350.513° 


7 Cru 


071 




266193414. 
183.104° 


266193267. 

183.241° 


266193260. 

183.248° 


266193238. 
183.267° 


266192889. 
183.578° 


266192521. 
183.886° 


7 Cru 


073 




267426088. 
182.136° 


267425942. 
182.281° 


267425935. 

182.288° 


267425913. 
182.309° 


267425565. 
182.640° 


267425200. 
182.969° 


7 Cru 


077 




269858216. 
181.086° 


269858071. 
181.240° 


269858064. 

181.247° 


269858043. 
181.269° 


269857699. 
181.619° 


269857336. 
181.965° 


7 Cru 


078 




270466692. 
180.860° 


270466548. 
181.015° 


270466541. 
181.023° 


270466520. 
181.045° 


270466176. 
181.398° 


270465814. 
181.747° 


/3 Gru 


078 












270512795. 
302.800° 


270512284. 
294.401° 


7 Cru 


079 




271045522. 
179.175° 


271045367. 
179.354° 


271045359. 
179.363° 


271045337. 
179.389° 


271044968. 
179.795° 


271044581. 
180.198° 


RS Cnc 


080 




271872473. 
90.240° 


271872088. 
87.858° 


271872070. 
87.749° 


271872018. 

87.431° 


271871265. 
82.984° 




RS Cnc 


080 


c 


271877226. 
121.515° 


271877611. 
123.897° 


271877629. 
124.005° 


271877681. 
124.324° 


271878434. 
128.770° 




7 Cru 


081 




272320388. 
178.322° 


272320233. 

178.510° 


272320226. 

178.519° 


272320203. 

178.546° 


272319835. 

178.974° 


272319448. 
179.397° 


7 Cru 


082 




272956171. 
177.862° 


272956016. 
178.056° 


272956009. 
178.065° 


272955986. 
178.093° 


272955617. 
178.532° 


272955229. 
178.967° 


RS Cnc 


085 






275057262. 
97.236° 


275057224. 
97.000° 


275057119. 
96.337° 


275055932. 
89.067° 


275055073. 

84.114° 


RS Cnc 


085 


c 




275059898. 
114.136° 


275059935. 
114.371° 


275060040. 
115.034° 


275061227. 
122.303° 


275062086. 
127.256° 


7 Cru 


086 


i 


275503697. 
176.829° 


275503542. 
177.033° 


275503535. 
177.043° 


275503512. 
177.073° 




275502756. 
177.995° 


RS Cnc 


087 


i 










276329999. 
92.710° 


276328974. 
86.701° 


RS Cnc 


087 


c 










276334140. 
118.549° 


276335165. 

124.557° 


7 Cru 


089 




277408751. 
176.576° 


277408596. 
176.781° 


277408589. 
176.791° 


277408566. 
176.821° 


277408199. 
177.287° 


277407813. 
177.749° 


7 Cru 


093 




280045204. 
208.249° 


280045028. 
208.061° 


280045020. 
208.052° 


280044994. 
208.024° 


280044576. 
207.598° 


280044136. 
207.175° 


7 Cru 


094 




280681410. 
191.683° 


280681250. 
191.696° 


280681242. 
191.697° 


280681218. 
191.699° 


280680836. 
191.728° 


280680433. 
191.758° 


7 Cru 


096 




282014259. 
185.190° 


282014112. 
185.280° 


282014104. 
185.285° 


282014083. 
185.298° 


282013732. 
185.504° 


282013362. 
185.708° 


7 Cru 


100 




285034037. 
224.282° 


285033857. 
223.835° 


285033848. 
223.814° 


285033822. 
223.749° 


285033398. 

222.741° 


285032956. 
221.750° 


7 Cru 


101 




285861190. 
224.289° 


285861011. 
223.842° 


285861002. 
223.820° 


285860976. 
223.755° 


285860552. 
222.746° 


285860110. 
221.755° 


7 Cru 


102 




286686360. 
223.942° 


286686182. 
223.500° 


286686173. 
223.479° 


286686147. 
223.415° 


286685725. 
222.418° 


286685285. 
221.438° 


13 Peg 


104 




288914432. 

342.574° 


288914336. 
343.021° 


288914332. 
343.042° 


288914318. 
343.107° 


288914091. 
344.113° 




RCas 


106 




291039691. 
90.705° 


291038969. 

86.728° 


291038939. 
86.566° 


291038853. 
86.097° 


291037730. 
80.186° 


291036830. 
75.723° 


a Sco 


115 




302022977. 
157.895° 


302022638. 
158.409° 


302022621. 
158.434° 


302022571. 
158.508° 


302021771. 
159.664° 


302020939. 
160.797° 
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TABLE 6 

Time, longitude and phase differences used to determine pattern speeds (full table included in on-line supplement) 



Wave 


Occultation Pair" 


5t (days) 


5X (dcg.) 


<5</> (dcg.) 


0-4, (dcg.) 


W82.00 


RSCnc085c-RSCnc085i 


0.03051 


16.9 


250.6 


3.0 


W82.06 


RSCnc085e-RSCnc085i 


0.03137 


17.4 


249.1 


5.7 


W82.21 


RSCnc085e-RSCnc085i 


0.03381 


18.7 


251.1 


7.7 


W84.64 


RSCnc087e-RSCnc087i 


0.04799 


25.9 


235.5 


4.2 


W80.98 


RSCnc080c-RSCnc080i 


0.05501 


31.3 


141.6 


5.0 


W84.64 


RSCnc085c-RSCnc085i 


0.06133 


33.3 


195.5 


8.9 


W82.00 


RSCnc080e-RSCnc080i 


0.06393 


36.0 


148.5 


8.8 


W82.06 


RSCnc080e-RSCnc080i 


0.06434 


36.3 


148.1 


6.1 


W82.21 


RSCnc080e-RSCnc080i 


0.06555 


36.9 


150.1 


4.3 


W87.19 


RSCnc087e-RSCnc087i 


0.07165 


37.9 


191.4 


12.1 


W84.64 


RSCnc080e-RSCnc080i 


0.08301 


45.8 


145.4 


8.5 


W84.64 


betGru078i-gamCru078i 


0.53957 


121.4 


47.3 


1.7 


W84.64 


gamCru081i-RSCnc080e 


5.10878 


50.2 


171.7 


7.2 


W82.21 


gamCru081i-RSCnc080c 


5.12178 


54.2 


208.5 


4.0 


W82.06 


gamCru081i-RSCnc080c 


5.12265 


54.5 


136.6 


5.9 


W82.00 


gamCru081i-RSCnc080e 


5.12294 


54.6 


110.6 


2.3 


W80.98 


gamCru081i-RSCnc080e 


5.12919 


56.8 


17.1 


11.3 



Each occultation is designated by the sequence star-name, rev number and ingress/egress. 



